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Bipolar spherical harmonics (BiPoSHs) provide a general formalism for quantifying departures 
in the cosmic microwave background (CMB) from statistical isotropy (SI) and from Gaussianity. 
However, prior work has focused only on BiPoSHs with even parity. Here we show that there 
is another set of BiPoSHs with odd parity, and we explore their cosmological applications. We 
describe systematic artifacts in a CMB map that could be sought by measurement of these odd- 
parity BiPoSH modes. These BiPoSH modes may also be produced cosmologically through lensing 
by gravitational waves (GWs), among other sources. We derive expressions for the BiPoSH modes 
induced by the weak lensing of both scalar and tensor perturbations. We then investigate the 
possibility of detecting parity-breaking physics, such as chiral GWs, by cross-correlating opposite 
parity BiPoSH modes with multipole moments of the CMB polarization. We find that the expected 
signal-to-noise of such a detection is modest. 



I. INTRODUCTION 



The detection of anisotropics in the cosmic microwave background (CMB) [T] has revolutionized the precision with 
which cosmological measurements can be made. Most of the information that has been obtained from the CMB so 
far has come from its power spectrum, the two-point correlation function, under the assumptions of isotropy and 
homogeneity. However, in recent years, attention has been paid to effects that go beyond the power spectrum, such as 
weak lensing [2], cosmic birefringence [31 0], and departures from statistical isotropy (SI) [IHZ] and from Gaussianity 

0DH). 

Bipolar spherical harmonics (BiPoSHs) [11H13) provide an elegant and general formalism for quantifying a number 
of these physical effects. If the CMB map is Gaussian and statistically isotropic, then its statistics are specified 
entirely in terms of the power spectrum Ci , the expectation value of the squared magnitude of the spherical- harmonic 
coefficients ai m for the map, and there are no correlations between different et/ TO s. A wide variety of departures from 
SI and Gaussianity induce correlations between different a; m s. The point of the BiPoSH formalism is to parametrize 
correlations between two different coefficients, a/ m and ai> m >, that represent two different "angular-momentum" states, 
in terms of total angular momenta L and M. Bipolar spherical harmonics have been used to search for non-standard 
cosmic topology [2] , anisotropy in primordial power [TS] , and model- independent departures from SI [131 [TSH^T] . 
They have also been used to test for asymmetric beams [22] and/or other systematic artifacts in WMAP [23 . BiPoSHs 
for polarization have been proposed to search for position-dependent rotation of the CMB polarization [2"ill26| . 

However, there is still more that can be done with bipolar spherical harmonics, and the purpose of this paper is to 
enumerate some of them. First and foremost, we point out here that almost all prior work on BiPoSHs has considered 
only BiPoSHs with even parity (Sec. [n]). There exists an entire other set of BiPoSHs that have the opposite parity, 
and these can provide probes of both cosmological effects and systematic artifacts that would remain elusive with the 
even-parity BiPoSHs that have been considered so far. We show, for example, that lensing by gravitational waves 
(GWs) can excite odd-parity BiPoSHs, and we describe a pointing error that could also excite these modes. In the 
process, we also show how gravitational lensing, by both density perturbations as well as GWs, can be described in 
terms of even- and odd-parity BiPoSHs (Sec. Ill I. Finally, we discuss how odd-parity BiPoSHs could be used as probes 
of parity violation, and consider in particular the cross-correlation of opposite parity CMB lensing and polarization 
components (Sec. IV). We calculate the anticipated spectra and errors for such correlations, and determine that a 
large signal-to-noise is not expected for these cross-correlations, given the current upper bounds on a GW background. 
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II. REVIEW OF BIPOLAR SPHERICAL HARMONICS 

A. Statistically Isotropic and Gaussian Maps 

A CMB temperature map Tin), as a function of position n on the sky, can be decomposed into spherical- harmonic 
coefficients 

ai m = J d 2 hT{h)Y lm {n). 

If the map is statistically isotropic and Gaussian, then the statistics can be determined entirely in terms of the power 
spectrum C;, defined by 

(ai m a*, m ,) — Ci 5u> 5 mm >, (1) 

where the angle brackets denote an average over all realizations, and 8w and S mm i are Kronecker deltas. Eq. 
states that all of the ai m are uncorrelated, and Gaussianity further dictates that the probability distribution function 
for any ai m to take on a particular value is a Gaussian distribution with variance C/. 

The spatial temperature autocorrelation function is defined to be C(n,n') = (T{n)T{n')) . Most generally it is a 
function of the two directions n and h' . However, if the map is statistically isotropic and Gaussian, then the spatial 
correlation function depends only on the angle 9, given by cos# = n-h', between the two directions. In this case, 

where Pi{x) are the Legendre polynomials. 

B. Departures from Gaussianity/SI 

Departures from Gaussianity and/or SI will induce correlations between different a; m s. The most general correlation 
between any two a; m s can be written, 

(<H m a1, m ,)=C l 5iv8 mm .+ ]T (-l) m '{lml', -m'\LM)Aff , (2) 

LM;L>0 

where Ci is the (isotropic) power spectrum, (I ml' m'\LM) are Clebsch-Gordan coefficients, and the A\ft l are BiPoSH 
coefficients. The spatial two-point correlation function is then 

C(n,n')=Y J ^^C l P l (n-n')+ ]T A*f {^(h) ® Y v {h')} LM , (3) 

l U'LM 

where 

{Yi{n)®Y v {n')} LM = ^{l 

ml' m'\LM)Yi m (n)Yi' m >(n / ), (4) 

mm' 

are the bipolar spherical harmonics (BipoSHs). These BiPoSHs constitute a complete orthonormal basis for functions 
of n and n' in terms of total-angular-momentum states labeled by quantum numbers L and M composed of angular- 
momentum states with Im and I'm'; they are an alternative to the outer product of the {I, m} and {I' , m'} bases. 

C. Odd-Parity Bipolar Spherical Harmonics 

It is instructive to decompose AfjJ^ into its odd and even parity parts, 

Af^=A® w 1 ^ -+A e w 1 J , (5) 
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where A®^ V M (y4 e ^, M ) are zero for the sum I + I' + L being odd (even). It follows from the symmetry C(n, h') = 



C(h',fi) that A^f™ (A e ^ vl ) are (anti) symmetric in I and We also infer that A® u ^ = (-l) M A®^, M and 



-M 



aQ lm: 

A IV 



= (-l) M+1 A e Z ;, . Thus, odd-parity BiPoSHs vanish for I = V . Prior literature has considered physical 



effects (e.g., nontrivial topologies [27] . SI violation [El [28]) that produce only j4®^, M , the even-parity BiPoSHs, and 



measurements have been carried out with WMAP data only for the A®\^ [16l ES] • In this paper, we consider also 



the odd-parity BiPoSHs A e „, 

Estimators for the BiPoSH coefficients (both the © and modes) can be constructed from a map of the CMB 
temperature field T(h), as follows: 

A™ = £ W^Wf 1 aZ P 4^ p (-ir'{lml', —m'\LM), (6) 

mm' 

and this estimator has a variance, under the null hypothesis (a SI Gaussian map), 

(a^a§^*^ = s LL , s MM , [s lT 6 VP + (-1)^'+%-,%] cr p cr p wr 2 w„\ (?) 

where a^ p = W\ a lm + a"™ se and (7™ ap = Wf C[ + N t are the temperature spherical-harmonic coefficients and 
power spectrum corrected for detector noise and finite resolution. The Gaussian detector window function, which 
encapsulates the effects of finite detector resolution, is given by Wi = exp [— Z 2 FWHM /(161n2)] , where 6>fwhm is the 
full width at half maximum of the detector. The instrumental noise contribution to the temperature power spectrum 
is given by 

Ni = 4^(NET) 2 



where / s k y is the fraction of the sky observed, NET is the noise equivalent temperature of the detector, and t bs is 
the length of time over which the CMB was observed by a particular survey. We notice that the variance in Eq. ^ 
vanishes for odd parity and I = I' = I = I', which is expected given that odd-parity BiPoSHs with I — I' vanish. 

The noise in any individual Ah, is large, and so a search for a statistically significant departure from zero in one 
or a handful of AhJ^s will probably not be too effective. It is better to consider specific models and/or parame- 
terizations for departures from SI/Gaussianity and then combine the AfjJ^s into a minimum-variance estimator for 
the SI/Gaussianity- violating parameters of those models. For example, Ref. [TTJ [T2] considered the bipolar power 
spectrum kl = Y^wm I^h /M | 2 as a parameterization for departures from SI. As another example, Ref. [15] combined 
Afj^s with L — 2 and I' = I, I ± 2 to derive minimum- variance estimators for the amplitude of an inflation-induced 
primordial-power quadrupole of the type considered in Ref. |30j . 



III. GRAVITATIONAL LENSING 
A. Gradient and Curl- Type Deflections 

Consider a statistically isotropic and homogeneous Gaussian temperature map T g (n) on the sphere, where h is a 
position on the sky. Now suppose that each point on the sky n has been deflected from an original direction h + A(n) 
so that the observed temperature is T(n) = T g {n + A) ~ T g (n) + A • V gT g (n). This deflection might come about 
cosmologically through weak gravitational lensing or may arise as an instrumental/measurement artifact (for example, 
if there are pointing errors). 

The most general deflection field A can be written in vector notation as 

A = V -0(n) + V -x n(h), (8) 

or in component notation, Aj = (V g) . (f>(n) + eij (V g) ri(n), 1 in terms of two scalar functions <f>(n) and Q(h) on the 
sphere, where V g is the angular covariant derivative on the unit sphere. In other words, the most general vector field 



Here, the Levi-Civita symbol on the unit sphere can be defined in terms of its three-dimensional equivalent as eij = — r^. The 
choice of sign here can be understood as the choice to have the spherical polar coordinates (8, <j>) form a right-handed coordinate system 
on the sky, since it will ensure that the basis vectors satisfy eg X e^, = 1. 
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on a two-sphere can be written as the gradient of some scalar field <j)(n) plus the curl of some other field f2(n). Weak 
gravitational lensing by density perturbations gives rise, at linear order in the lensing potential, only to the gradient 
component. A curl component can arise cosmologically from second-order terms in the deflection field or from lensing 
by GWs. Systematic measurement effects may conceivably give rise to both types of deflections. 

We now show that the and A e ^/ M BiPoSH coefficients are induced, respectively, by the gradient and curl 

components of the deflection field. The change in the temperature moments induced by lensing is (at first order in 4> 
and O), 



<)(i, ni = I d 2 hY* m (n) | 



V e -< 



VgT{h) + V e -0(n) 



-T(n)]} 



~- E E a '' 

LM;L>0 I'm' 



where in the second line we have decomposed 



LM 



^gY LM {fi) ■ VgYi, m ,(h) 



LM 



^gY LM (n) 



in)] } 



l) = E Y LM(n) <j>LM, 

L=l M=-L 



(9) 



and similarly for f2(n). We do not consider L = modes of <j> and f2 since they would not cause a deflection. In the 
notation of Ref. |3T] , 



lm 



[{Yim m+ - -{Y lm m_ 



where {Yi m and _iYj m are spin-weighted spherical harmonics, the null coordinates rh± — (eg + ie^)/^/2, and the only 
non-trivial products of the null coordinates are m + • m_ = 1, and rh + x m_ = i. Thus, it is obtained that 



VgY LM ) x (VgY Vm > \ = 



\J L{L 4 
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W L ( L - 


f 1)Z'(Z' + 1) 



[{{Ylm) i-iYi'm') + {-iY LM ) (iYi> m i)\ , 

[(iYlm) {-lYv m >) - {-lY LM ) {lYvm,)] . 

Using the triple integral |31j of spin-weighted spherical harmonics, the 8ai m for the gradient and curl terms are 
obtained as 



Sain 



where 



E E 



LM:L>0 I'm 



(-l) M + m a llm , Gfj, 
sJ(2L +1)1(1 + 1) 



'LM ~ 



1 + (_ 1 y+l'+L 



in 



LM " 



i - (-ly+i'+L 



(l ml', —m'\LM), 



g\ v , ^ { L + m + mv^mi + mv + i ){ml l|Ll) _ 

Up to linear order in the deflection coefficients 4>lm and SIlMi the even- and odd-parity BiPoSH coefficients arc then, 



A® LM 



A e LM _ 

Ji ui — 



tlm 



V2L + 1 

iQlm 
V2L + 1 



y/V{V + 1) VW+V 

/~i r*L s~i r^L 

^l+'l'l ^l ,Lr W 



VW' + 1) VW+V 



— Qfi^LM, 



(10) 
(11) 



where we have defined the quantities 



i 



QL 



V2L + 1 



V2L4- 1 



CiG 



C V G 



y/VJjf + 1) v^^l + l) 

^}}fin ^i'^w 

VW' + 1) y^TT) 
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Clearly, the gradient part contributes only to and the curl part only to A e a , . Further, it is explicit that the 

gradient and curl parts of the deflection correspond, respectively, to the symmetric and antisymmetric (in {11'}) parts 
of the total A^f 1 . 

Suppose the A^s have been measured using the estimators in Eq. (|6|). If we then assume that lensing is the 



dominant source of BiPoSHs we can use Eqs. ([7]), ( [XO] ), and ( 11 ) to construct maximum- likelihood estimators for the 
gradient and curl components of the deflection field, 

0LM = t , (l^j 

y,u>\q?i-\ 2 / (wr'w^crcrn 

The variance of these estimators, under the null hypothesis of no lensing, is given by 



4>L>M> ) = 5lL> $MM> (^l) 2 = 2 $MM 



(14) 



EI^ L r/(^r 2 ^ 2 ^ map cr 



(15) 



where the sums in Eqs. (12) and (14 1 only include pairs of I' for which I + I' + L is even, while those in Eqs. (13) 
and (15) only include pairs for which this quantity is odd. 



B. Deflection Field from Metric Perturbations 



Cosmic shear, weak gravitational lensing due to density perturbations or GWs along the line of sight to the CMB, 
will produce displacements like those in Eq. (|8| . Our goal here will be to calculate the displacement spherical- harmonic 
coefficients 4>lm and Qlm that arise from gravitational lensing due to density perturbations and GWs. There is a vast 
literature on lensing by density perturbations and also specifically on lensing of the CMB by density perturbations 
[2] ■ Our density-perturbation results follow most closely those of Refs. [3U E2] ■ Lensing by GWs has been considered 
in Ref. [33]. We follow primarily the approach of Refs. |3_U[35j, who calculated Qlm due to GWs, but extend their 
results to include §lm from GWs, reproducing the results of Ref. [35]. We make use in this Section of relevant work 
on lensing and/or differential analysis on the celestial sphere in Refs. [3T1 1321 1371 135] . 

We write the metric for the perturbed spacetime as 



ds 2 — a 2 (rf) [—dr] 2 + (8ij + hij) dx*dx^ , 

where hij is the metric perturbation in the synchronous gauge, and r\ is the conformal time. Now consider a photon 
that we observe to come from the direction n on the sky. In the absence of perturbations, this photon travels along a 
path x(rj) = (?7o — rj) n as a function of conformal time t], where 770 is the conformal time today. Metric perturbations 
will induce perturbations in this trajectory, which we can calculate by integrating the geodesic equation back over 
the photon path to find the direction of propagation of the photon when it was emitted at a conformal time 77. To 
first order in the metric perturbation h, we find the original direction of propagation of the photon on the sky to be 
h + A, where [35] 

(16) 

[»?', (»)o-»)')™] 

Here, we have ignored the observer terms h^^io), and we have defined the projection tensor P im — <5 jm — n^rim onto 
the space perpendicular to the unit vector n. The subscript indicates that the quantities in the integral are evaluated 



A'(n) 



p. rv 
-=^- / d V ' 



h m jhj - -{rf - 7])n k fiid m hki 



() 



at time and space coordinates (rj, x) — [rf , (i]q — rf)n); i.e. they are evaluated along the unperturbed path of the 
photon. In our case, the source is the CMB, and 77 = ?7i ss is the conformal time at the surface of last scatter. However, 
the calculation could also be applied to the lensing of galaxies in which case the relevant conformal time would be 
that corresponding to redshifts z ~ 1. 

The functions (f){n) and £l(n) in the decomposition in Eq. ([8]) can be obtained from 



VU(h) = V r A(n), Vlfi(n) = -V* x A(n), 



(17) 



where as before V g is the angular covariant derivative on the unit sphere. As Ref. [25] notes, the standard lensing 
convergence is n = — (1/2)V| <j) and the lensing rotation is to — (1/2)V- 57. 
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The gradient component is obtained from 



V 2 s 4>(n) =V r A = - ; 



'/o 



dv'ivo ~ v'){$ik - hih k ) 



d k (hijn 3 ) +-(?/- rj)didk{hi m hih m ) 



Vo-V [J n 

+ I drj' [Zhihjhij - ha + (rj' - v) (^jdihij - 2 hihjfikdkhij)]^, ^^,^ >, 



In', (vo-v')fii 

(18) 



where we have used the fact that Vg, which acts on the unit vector n, behaves as V~ = (770 — f]')(^ik — fiink)dk inside 
the integral due to the dependence of x on n as defined in the integrand subscript. 

Let us now consider the curl component. For this calculation we must use V?.f2 = -V? x A and then note that, as 



before, V^. = (?yo — w')($ik — hihk)dk inside the integrand. Applying this to Eq. (16), we have [34 



VgO,(h) = - / drf (riimeijkdjhki)^, fi^o-n')]- 



(19) 



C. Lensing by Density (Scalar Metric) Perturbations 



Let us first consider scalar perturbations. In the conformal-Newtonian gauge in the absence of anisotropic stresses, 
the metric is given by 



ds 2 = a 2 (r]) [-(1 - 2$)d?7 2 + (1 + 2$)6 ij dx i dx i ] . 

Noting that a conformal transformation preserves null geodesies, our calculations of the photon path will be unaffected 
if we work in a synchronous metric obtained from the conformal-Newtonian form through multiplication by (1 + 2$). 
Assuming that $ is small and keeping terms only to linear order, we find the conformally related metric, 



ds 2 



a 2 (77) [-drf + (1 + 4$)%cfa;W] 



Using this metric perturbation hij = 4<!><5j J in Eq. ( 18 1 above, we find that the first, third, and fourth terms vanish, 



giving for the gradient-type lensing caused by scalar perturbations, 



2 f Vo 

V| <f ca (n) = / drftf - 77) [ (§ij ~ Mi) (Vo - r)')didj® -2^8^]. 

VO - V J q 

For small-scale fluctuations, the second term will be negligible compared with the first, so it can be dropped. We can 
rewrite the spatial derivatives in terms of to find 

V|0 sca (n) = — r dr fHLz^^ 9 (rf ) (no-^n), 

9 Vo-V J v Vo-V 9 

and we can remove the angular derivatives to obtain the usual expression for the projected potential 
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rf-n 



^(h) = -2 / drf- -<&(,/, ( Vo - V ')n). 

Jr, KVo ~ V){Vo ~ V ) 

We can once again decompose 4>(n) in terms of its spherical-harmonic coefficients as in Eq. We then find 

4>T M = J d 2 nYl M (n)^(n) 

= -2 r drf- tp- ^ [ d 2 h Yl M {n) $(r/, % - rf)n). (20) 

Jr, KVo -V){Vo -W) J 

Thus, lensing by density perturbations w ith a given projected potential is characterized by nonzero even bipolar 



spherical harmonics A®^, M given by Eq. (10 1 with 4>lm given by (jf^fj above. Scalar perturbations cause no curl- 



type len sing, which we can see in several ways. For scalar perturbations, hij oc $ 5y, and so the left-hand side of 
Eq. (19 1 vanishes. Then, by taking a Laplacian of the mode expansion Qlm = / d 2 n Q(n)Y£ M (n), and noting that 



the spherical harmonics are eigenfunctions of the Laplacian with eigenvalue L(L + 1), we can write 



Qlm = L(L + 1) J d ^ Y LM(n)Vp(h). (21) 

Thus, we find that all of the £I s £mi exce Pt possibly for the unphysical L — mode, vanish. Equivalently, an argument 
can be made that scalar perturbations have no preferred direction, and so could not generate curl-modes, which do 
have a preferred direction. Thus, scalar modes produce no odd bipolar spherical harmonics A e ^/ M . 



We can go on to find the autocorrelation power spectrum of the Starting from Eq. (20), we use the fact that 

the potential perturbations $(77, k) today are related to their primordial values $p(k) by 

Hv,k) = ^p(k)T^(k)^, 

where a(r}) is the scale factor, T sca (fc) is the scalar transfer function that describes the evolution of scalar modes 
through the epochs of horizon crossing and matter-radiation equality, and D\ (77) is the growth function that captures 
the scale- independent evolution of scalar modes at later times [4"0] . The transfer function can be approximated using 
the fitting form of Ref. [H] , 

T sca (a; = k/k cq ) = ln(1 Q ^ 7 ° 1 1 J a;) [1 + 0.284 s + (1.18 a;) 2 + (0.399 x) 3 + (0.490 a;) 4 ] ~ a25 , 

where k cq is the wavcnumbcr of the mode that crossed the horizon at matter-radiation equality, defined as 

k cq = a cq H(a cq ) = y/2 Hq a cq ^ 2 . We can write the growth function, under the assumption of cosmological- 
constant dark energy, as 

_ , . bn m H( V ) f a ^ da' 



3 • 



2 Ho J (a'H(a')/H ) 

We also write the autocorrelation of the primordial scalar fluctuations ($p(fc) $p(fc')) = (27r) 3 5 3 {k — fc')P$(fc), where 
the primordial power spectrum is given by 

50tt 2 / k \™ 3_1 a2 / ft 



9k 3 \H J H \Di(a = 1) 

With these ingredients, and after using the partial-wave decomposition, 



e ik{i}a-v') < 



= J2 * L ( 2L + ( k (m - rf)) Pl(cos0), (22) 



L=0 
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FIG. 1: Here we plot the autocorrelation power spectrum C"£ of the gradient-type <f> modes of cosmic shear. In green squares 
we show the autocorrelation of the <f> modes from lensing by scalar perturbations, and in blue circles that of the <j> modes induced 
by tensor perturbations. We use the WMAP-7 cosmological parameters, and assume the maximum allowable tensor-to-scalar 
ratio r = 0.24 from the WMAP-7 data combined with BAO and the Ho measurement [42], to calculate the tensor contribution. 
The error with which these power spectra could be measured using the parameters of the Planck satellite is shown as red +s. 



we find the autocorrelation power spectrum to be 



CI 



5 (?7 - m) 



dk k 2 P^(k)T sca (kf 



(rjo - rf) a(r)') 



(23) 



To calculate the magnitude and shape of this autocorrelation function, we employ the WMAP 7-year cosmological 
parameters of Ref. [32] . We plot the result of our calculation in green squares in Fig. [l] 



D. Lensing by GWs (Tensor Metric Perturbations) 

If the metric perturbation hij is caused by GWs, we can decompose it into plane waves, 

/,/3l _ 
_ e *-*T(M) E fc°W ( 24 ) 
^ n > a=+,x 

where we sum over the two GW polarizations + and x, the plane- wave amplitudes are h a (k), and eg- are the polar- 
ization tensors, which are transverse, traceless matrices. Here, T(k,rj) is the GW transfer function, which gives the 
conformal-timc evolution of the mode; Ref. |34j notes that it is well approximated by T(k,rj) = 3ji(kr))/(kr)). 

Now consider a single GW mode propagating in the z direction with wavenumber k and + polarization. In this 
case the polarization tensor is 
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FIG. 2: Here we plot the autocorrelation power spectrum C L of the curl-type modes of the weak lensing of the CMB 
temperature field. These modes can only be induced by tensor perturbations. We show the signal in blue circles and the error 
with which they could be measured using the parameters of the Planck satellite as red +s. 




efj(kz) = 

The only non-zero metric-perturbation components are then h xx = —h yy = h + (k) e tkz T(k,n). The unit vector 
n = (sin 9 cos <p, sin 9 sin ip, cos 9). The curl component of lensing of the CMB by tensor perturbations is then 

\7jn tcn (h) =ikh+(k) sin 2 9 sin 2(p / drf T{k, j/Je^*-" ) c ° s «. (25) 

J 7} 

A GW with the x polarization is the same as that with the + polarization, but rotated by 45° to the right. The 
f2 ton (n) pattern is therefore the same, but with sin2<p replaced by — cos2</9. We thus see that lensing by GWs will 
give rise to nonvanishing A Q f^ . 

The gradient component of cosmic shear due to tensor perturbations is a bit more complicated; it is 



^ sin 2 9 cos 2(p J drf T(k, rf) |3 - 2ik{rf - 77) cos 9 



v 2 r tcn (n) = -— {k) 
ti no 

+ (Vo - V') 

This can be further simplified by noting that 



k 2 

ik cos 9 — —(r]' — rj) sin 2 9 



e ifc(j?o-»/) cose 



ik cos <9e 4fc(, '°-'>' )cose = JL e iHvo-v')<y°*(> 
on' 
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which then leads to 



Vf</> tc » 



h + (k\ r no ( 8 

^ sin 2 9 cos 2<p / drf T(k, rf) \ 3 + 2 (W - rf) 



(vo - rf) 



drf 
d 



i W y) ( 2 a 2 

dr/' 2 V dr/ 2 



a ik(ri -7j') cost 



(27) 



For the x polarization, we replace cos 2^5 by sin 2ip. 

Note that the expressions for V 2 (j) ten and V 2 f2 tcn differ only in two ways: (1) The curl mode has a sin2<y9 
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dependence on the azimuthal angle ip, while the scalar mode has a cos2<p dependence (for the + polarization). (2) 
The rj' dependences of the two integrands differ. 

We now find the spherical-harmonic coefficients </>^ = / d 2 n (/> ten (n)Y£ M (n) and 0^ = J d 2 hft ten (h)Yl M (n). 



Taking the angular derivatives of this decomposition of the curl component, we find the result Eq. (21). We also 
expand these coefficients in terms of their polarization and k modes, 



r-)ten 

"LM — 



d 3 k 



(2tt) 3 

y ' a=+ 



(28) 



If we consider just one mode, with a 
a conformal time integral: 



and k = kz, and use Eq. (251, its amplitude simplifies into an angular and 



ng° + (fcl) = - tk ^ k \ r° drf T(k,r]') [ d 2 nY£ M (n)sm 2 sin(20) e iHvo-v')co S e 



L(L + 1) J n 

The azimuthal integral is easily taken once the spherical harmonic is decomposed as 



Y L\i(ri) 



'2L + 1 (L - M)\ ^_ 
4?r (L + M)\ " 



iM<t> 



P LM (cos ( 



and yields the result that only M = ±2 modes remain. The polar integral can then be taken by using the partial-wave 
decomposition Eq. ( 22 ) and by converting associated Legendre polynomials into regular Legendre polynomials and 



using their orthogonality. The final result that we obtain for the spherical-harmonic coefficients of the curl mode is 



Qtcn - 
5 LM 



(kz) = i L h+{k) (5 Mi2 - Sm.^J 21 ^ Ff(k), 



where 



F£(k) 



l2n{L + 2) 



1 



rk,1 ° , rpi J L (kvo-w) 
, dw 1 (w)- 
(L-2)\ L(L + l)J kv K 



(fcryo - w) 



(29) 



(30) 



is a transfer function for fL Note that in writing Eq. ( [30] ) we have assumed that T(k,rj) = T(krj), and that for the x 
polarization the sin 2^? dependence of Q(n) is replaced by — cos2ip, so that the factor (5m,2 — is replaced by 

-i(5 M ,2 + §M,-2)' 



Likewise, noting the similarities between Eqs. (25) and (27), and decomposing 



'LM 



into modes as in Eq. ( 28 1 



d z k 
(2tt) 3 



(31) 



Ct = + ,X 

the result for the amplitude of the gradient mode with a = + and k = kz is 



0LM + (kz) = l L h + {k){8 M ,2 + S M ,- 2 ) 



2L + 1 



Ft(k), 



(32) 
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where 



I2it(L + 2)\ 



(L-2)\ L(L + l) 



dw u — 



w 



V) 



T(w) 



d^ + 2 {w - k7l) 



d 2 
dw 2 



jL(ki] Q - w) 



w 



(33) 



Again, the factor (Sm,2 + <5a/,-2) is replaced by — i(5m,2 — &m,-2) for the x polarization. 

The contributions from this Fourier mode to the 4> tcn and fi ten power spectra are Cf 



\kz) 



Em<I^lmI 2 )/( 2L + !) and c2 ntcn {kz)+ = EA/<l fit LMl 2 ) /( 2L + !)■ Note that it; is onl y the M = ± 2 modes 

that contribute. By rotational invariance, the contribution from the x polarization is the same, as is the contribution 
from any other mode with the same wavenumber k but pointing in a different direction. If the gravitational waves 
have power spectrum Pr(k), defined by 



(2TT) 3 S D {k-k')S lJ P T {k), 



(with = {x, +}), then the <p and ^ power spectra are 

gxxtcn ■ i d 3 k 



(2tt)' 



P T (k) [F?(k)] 



(34) 



(35) 



for X = {0, In this paper, we will assume a scale-invariant power spectrum 



P T (k) 



2fc 3 



A 



(36) 



where we have neglected the spectral tilt and adopt the parameters of WMAP7 . 

We calculate the variance in the measurement of these autocorrelation functions from an observed CMB temperature 
map, under the null hypothesis of no GWs, and obtain an expression in terms of the variance of the <p and ^ estimators, 
Eqs. (fl4l) and tfl5] 



AC; 



2L 



C<p<p sea 
L 



(37) 



ACf = 



.02 



2L + 1 



(38) 



Here, cr^ 2 and cr^ 2 are the variances of our estimators for <\>lm and Q,lm as found in Eqs. (14) and (151. Note that 



under the null hypothesis of no GWs, there is no expected cosmological curl-type lensing signaT7so this term is absent 
in Eq. ((381). 



To calculate these autocorrelation functions and their variances, we use the WMAP 7-year cosmological parameters 
of Ref. [32]. We plot the results of our calculations in Figs, [l] and [2j where lensing from scalar perturbations is 
plotted in green squares (absent for since there are no scalar contributions to the curl modes), lensing from 

tensor perturbations is plotted in blue circles, and the variance of these measurements is shown in red +s. We use the 
parameters of the Planck satellite, NET = 62 fiK s 1 / 2 , t ohs = 2yr, 6> FW hm = 2.0635 x 10~ 3 rad, and / s k y ~ 1- We can 
see that the scalar <f> signal is several orders of magnitude greater than the tensor signal, and that the variance with 
which the 0-</> power spectrum could be measured with Planck is higher than the scalar signal for low multipoles. The 
corresponding variance with which the power spectrum could be measured is also significantly larger than the 
signal. In both of these cases, therefore, the signal to noise of measuring the lensing from tensors using all multipoles 
with Planck is negligible, and remains negligible even in the case of the ideal CMB experiment with zero noise. 

Thus, a stochastic background of GWs with power spectrum Pr(k) predicts a spectrum of nonzero even and odd 
BiPoSHs given by Eqs. (10) and (11), with values of 4>lm and Qlm selected from Gaussian distributions with the 

r~<d><b ten 

variances GT 



and C2 ntcn 



given by Eq. ( 35 ) 



E. BiPoSHs from Pointing Errors 



A telescope pointing error can be described as a process that causes the positions of points on the sky to be 
mislabeled. This then causes an effective deflection of the points on the sky n Q b s = n + A(n), where n Q b s is the 
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direction that the telescope believes it is pointed in and n is its actual pointing direction. As we saw in Sec. Ill A 



we can decompose this deflection field A (ft) into gradient and curl components, which source even- and odd-parity 
BiPoSHs, respectively. Thus, from Eq. ( [l7| we can see that any pointing error that has a nonzero curl component 
Vg x A(n) will excite odd-parity BiPoSHs. 

Imagine, for example, that a satellite such as Planck misestimates the rate with which it is precessing. Since it is 
this precession that builds up observations of subsequent rings of the sky, such a misestimation would cause a shearing 
of each ring relative to its neighbors. This type of a deflection has a nonzero curl component, and thus would excite 
odd-parity BiPoSHs. Measurement of these BiPoSHs, and in particular the odd-parity BiPoSHs, can therefore provide 
a useful check for such pointing errors. 



IV. BIPOSHS AS PROBES OF PARITY VIOLATION 



A. Correlation of Opposite-Parity Lensing Components 



Since the A® w and A e ^f have opposite parity for the same L and M, a cross-correlation between the two can arise 
only if there is some parity-breaking in the physics responsible for producing the departures from SI/Gaussianity. Here 
we mention, by way of example, chiral GWs as a mechanism to produce such a parity- violating correlation |43H45| . 

The contribution to the cross-correlation power spectrum from a single Fourier mode in the z direction with + 
polarization is (kz) — Yl m ^lm) /(^L + 1) = 0; it vanishes as the contribution from M = 2 is canceled by 

that from M — —2. And if this is true, then by rotational invariance it is true for any other linearly-polarized GW. 
We thus conclude that a stochastic GW background predicts Cf^ = 0. In other words, there is no cross-correlation 

between <j> and f2, and thus no cross-correlation between the even and odd BiPoSHs, A®^ 1 and A Q f^ . 

Following Ref. [44] . however, consider a right-circularly polarized GW: hp = h + +ih x (i.e., we sum a + polarization 
wave with a x polarization wave out of phase by 90°). The azimuthal- angle dependence for the wave is then e 21 ^, 
and VLlm and 4>lm have contributions only from M = 2. There is thus a nonzero cross-correlation between <f> and 0. 
Similarly for a left-circularly polarized GW = h + — ih x , the ip dependence is e _2 * v , and only M = —2 modes are 
excited. There is again a cross-correlation between <fi and J7, but this time with the opposite sign. 

In the standard inflationary scenario, there are equal numbers of right- and left-circularly polarized GWs, and the 
cross-correlation between <j) and £1 therefore vanishes. But if for some reason there is an asymmetry between the 
number of right- and left-circularly polarized GWs 43-46 , a manifestation of parity breaking, then there may be a 



5 LM 



parity- violating cross-correlation between <f> and O, and thus between A®^ 1 and . 

The chirality of the GW background can be parametrized by an amplitude A which can take values between — 1 
and 1, where A = +1 denotes that all of the GWs are right-circularly polarized, and A = — 1 denotes that they are 
all left-circularly polarized. But we have seen that a right-handed GW contributes only to M = 2 modes, while a 
left-handed one contributes only to M = —2. We can denote this by weighting M = 2 components by (A + l)/2 and 
M = —2 components by (A — l)/2, so that our version of Eq. (29), for example, that is appropriate to the case of a 
chiral GW background will be 



^ 0M.2 ~ Ojw 2 



2L + 1 



(39) 



and similarly for Eq. ( 32 ) . In this way a fully right circularly-polarized GW background will have only contributions 
from M = 2, a fully left-circularly polarized background will have only contributions from M = —2, and if the amount 
of left and right-circularly polarized waves is equal, that is if the GW background is non-chiral, the contributions from 
M = 2 and M = —2 cancel. The 0-f2 cross-correlation power spectrum is given by 



C 



A 



d 3 k 
J2wf 



P T (k)Ft(k)FS(k). 



(40) 



Refs. [34, 36, 47] have shown that the amplitude of the stochastic gravitational-wave background is probably too 
small, even with the most optimistic assumptions, to produce a detectable gravitational-lensing signal in the CMB. 
The example of a chiral gravitational-wave background as a possible source of a detectable parity-breaking BiPoSH 
correlation is principally of academic interest. Still, Ref. [35] has recently argued that weak lensing of the CMB by GWs 
may be detectable in its cross-correlation with the CMB-polarization pattern induced by these GWs [37j [38J I48H50] . 
We thus surmise that a chiral gravitational-wave background may still be able produce a detectable parity-breaking 
signal in BiPoSHs in cross-correlation with the CMB polarization, an idea we explore in the next section. 
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B. Large-Angle CMB Polarization Spectra 

We follow the work of Ref. [35] , finding the multipole moments of the CMB E- and B-type polarization spectra for 
large angular scales by considering only those modes that are produced after reionization. The spherical-harmonic 
coefficients of B-type polarization modes can be decomposed as 

B ^ = / 72^3 E BM)> ( 41 ) 

where Bf m (k) is the amplitude of polarization B modes multipole moment Im in the direction k. The general form of 
this amplitude is quite complicated, but we can simplify it if, as in Sec. |Iff D| we consider only a single, +- polarized 
GW traveling in the z direction with wavenumber k. In this case, the B-mode amplitude can be written 



B+Jkz) = i l h a (k) (5 m , 2 - 6 m ,- 2 ) \[^F l B {k) 1 



(42) 



1 



= ^ [°dvm {{l + 2)j«-i[Afa, - f))] -(I- l)3l+i[Kvo v)}} fj dx 3 x y , 



(43) 



scattering rate 
and a the scale 



where the h a (k) are the amplitudes of GW modes as defined in Eq. (24), f{rj) is the 
f(rj) — n e (j]) gt o,(r]), with n e the electron density, ut the Thompson scattering cross-section, 
factor, and r) le and ry are the conformal times at reionization and today, respectively. Since we are only interested in 
small scales, we find the approximation r/i ss — is sufficient for our purposes, making the last integral significantly 
faster to evaluate. The result above agrees with the results of Ref. [35], whose method we followed in its derivation, 
up to a factor of i. 

We find that the corresponding E-type polarization multipoles from tensor perturbations take the same form as 
Bi m above, except for the opposite sign in front of S mt - 2 and a different factor in the curly brackets in Eq. M3J). 
From Ref. [49] we find this alternative form to be (21 + l)/2 {— ji(x) + j['(x) + 2ji(x)/x 2 + Aj[(x)/x\, where here 
x = [k(rjo — T))], and derivatives are with respect to x. Employing spherical Bessel function identities, we can then 
write the E-type polarization multipoles as 



^=/(0 q E x ^)' ( 44 ) 



E+Jkz) =z l h a (k) (6 m , 2 + «5 m ,_ 2 ) ^^±i F t E (k), (45) 
Z(Z + 3) . n . (l + l)(l-2) . n . ..1 f kr > -3j 2 (x)j 2 (kr)-x) 



' 2(2/-l) JMl wu In 2(21 + 3) ">$J kmss x (kri-x)* 

where the terms are defined as they were above for the B\ m amplitudes. 



C. Parity- Violating Correlations from Chiral GWs 



We now want to calculate the expected cross-correlation between CMB-polarization multipole coefficients and weak- 
lensing-induced BiPoSHs of opposite parity. Note that these cross-correlations are directly related to the parity-odd 
three-point correlations discussed in Ref. [51) . As we mentioned above, if there is no parity- violating physics, then 
in the cross-correlation of a parity-even and a parity-odd observable, M — 2 terms and M — —2 terms will cancel 
each other, giving a net zero cross-correlation. However, if for example the GW background is chiral, then parity is 
broken and we can get a non-zero cross-correlation between opposite parity observables. As we saw in Sec. |1V A[ a 
right-handed GW contributes only to M = 2 modes, while a left-handed one contributes only to M = —2. If we carry 
out a similar procedure for Eqs. (42), and (45) as we did in Eq. (39), weighting M = 2 components by (A + l)/2 and 
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FIG. 3: Here we plot the cross-correlation between the gradient (f> modes of the weak lensing of cosmic shear with the 
curl-type B modes of the CMB polarization in blue circles, and the noise on this measurement due to cosmic variance and 
Planck satellite instrumental noise in red +s. Since these quantities are of opposite parity, in the absence of parity-breaking 
physics we expect this cross-correlation to vanish. However, if we assume for example that the entire allowable GW background 
is right-circularly polarized, such a cross-correlation could occur. The cross-correlation is linearly proportional to the chirality 
parameter A, defined such that A = 1 denotes a completely right-circularly polarized GW background, A = —1 denotes 
completely left-circularly polarized, and A = denotes an unpolarized background. Here we assume the maximum allowable 
tensor-to-scalar ratio r = 0.24, the limit from WMAP-7 data combined with BAO and the Ho measurement |42j . Cusps in the 
absolute value of the correlation function correspond to sign changes of the correlation function. 

AI = —2 components by (A — l)/2, we can calculate parity-violating correlations between polarization and lensing 
components while accounting for the amplitude and handedness of a chiral GW background. 

First considering the cross-correlation between S-modes of the CMB polarization and gradient-type modes of cosmic 
shear, we write 



2L 



M 



As before, by rotational invariance we know that both + and x polarizations will contribute equally to C^ S , as will 
modes with any wavenumber k whose magnitude k is the same. We can see that only (^lm wm contribute to this 
correlation, and not 4> s t c m, as the scalar perturbation field is not correlated, on average, with the tensor perturbation 
field. Then using Eqs. (31), (32), (34), (41) and (42), we can write this cross-correlation as 



A 



d 3 k 
(2tt) 3 



P T {k)Ft{k)FE(k). 



(47) 



Similarly, we can write the cross-correlation between _E-type polarization modes and curl-type modes of cosmic 
shear, using Eqs. (p8l, (|29|, (1341), |44l, and (|45|, as 
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FIG. 4: Here we plot the cross-correlation C L between the curl-type Q modes of cosmic shear with the gradient-type E- 
modes of the CMB polarization in blue circles, and the noise on this measurement due to cosmic variance and Planck satellite 
instrumental noise in red +s. As with the cj> — B correlation, we assume a completely right-circularly polarized GW background, 
with the maximum currently permitted tensor-to-scalar ratio. 



CS E = A I ^3 P T (k) F?(k) FE (k), (48) 



where the GW power spectrum is given by 



2fc 3 

We want to calculate the magnitude and shape of such correlations, to determine whether such a signal is observable. 
We use the WMAP 7-year cosmological parameters and assume the maximum allowable level of GWs from early 
universe physics, with a tensor-to-scalar ratio r = 0.24, the limit from the WMAP-7 data combined with BAO and 
the Ho measurement [32] ■ We also assume that the GW background is entirely right-circularly polarized. As a first 
estimate, we calculate the level of such correlations while making several assumptions. We use the approximate form 
of the GW transfer function T(k, ry) ~ 3ji(krj)/(kr)), assume that reionization happened instantaneously so that the 
electron density n e is equal to a step function, and neglect contributions to the polarization modes that came from 
last scattering. The two last assumptions affect mostly the higher-!/ multipoles, which in this cross-correlation are 
suppressed since we see that ^lm anc ^ ^Im f &u °ff ver Y f as t with L. 

With these assumptions, we have calculated the correlation functions Cf 3 and C£ B , and show them as the blue 
circles in Figs. [3] and [4] Note that the absolute value of the correlation functions are plotted, and that the cusps in the 
profiles result from sign changes. Note also that both correlation functions are linearly proportional to the chirality 
parameter A, so that they would flip in sign if the GW background were left instead of right-circularly polarized. We 
are only interested in low multipoles, since our assumptions break down for larger L, and such multipoles are strongly 
suppressed in correlation with the weak-lensing modes. 
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D. Variance of 0-B and f2-E Correlations 

It is useful to know the variance with which we could measure such parity-violating cross-correlations. From Rcf. [37 
e see that the variance with 
variables X and Y is given by 



we see that the variance with which we could measure the cross-correlation C XY of two distinct Gaussian random 



where C XY — 1/(2L+ 1) Y2m XlmY£,m is the estimator for the cross-correlation, and C XY is its theoretical value 
under the null hypothesis. Ref. (37) then evaluates this variance, assuming distinct X and Y , to be 



1 



(AC? Y Y = ^ (C XY f + C xx ™* Cl Y ™ P , (49) 



IIC 



and N xx the noise 



where, as before, C XX map = W^Cl + N xx , with Wl the window function defined in Sec. 
in the measurement of C xx . 

In our case, the null hypothesis is that there is a GW background with the maximal tensor-to-scalar ratio, but it 
contains equal numbers of right- and left-circularly polarized GWs, i.e., it is not chiral. In this case, the theoretical 



value of parity- violating cross-correlations is zero, so that the first term in Eq. (49) vanishes. Then, assuming that 

4>lm an d Qlm are Gaussian random variables, a reasonable assumption since many uncorrelated noise processes are 
likely to contribute to this measured value, we find for the variances, 

1 



(AG*")' 
( 



2L + 1 



£r</)</)map ^jBBra&p (50) 



^£if2is\ 2 \ pfiS! map qEE map (51) 



To calculate these errors, we know that the instrumental errors on the polarization power spectra are given by 

EE _ AT BB _ 87r(NET) 2 



NZ» = N L 



We use the Planck-satellite parameters, as in Sec. HID We also use the CMB anisotropy calculator CAMB to calculate 
the temperature and polarization power spectra including effects at the surface of last scatter |52j . The resulting errors 
are shown as red +S in Figs. [3] and [4] This noise, which combines instrumental and cosmic- variance sources, is at 
least an order of magnitude above the corresponding maximum signal level at low multipolcs, and drops less rapidly 
with I so that the low multipoles yield the highest signal-to-noise. 

E. Signal-to-Noise Ratio of Chiral GW Background Detection 

We finally wish to calculate the achievable signal-to-noise of a measurement of the magnitude of such cross- 
correlations given our calculations of their shapes and variances. Such a measurement would tell us about the 
presence or absence of a chiral GW background, or of parity violation in the processes that caused departures from 
Gaussianity/SI in general. We can phrase the aim of this calculation as finding the error with which we could measure 
the chirality parameter A, which sets the amplitude of the cross-correlations relative to their maximum values in the 



case of a completely circularly polarized GW background, as in Eqs. (47) and (48). Let us calculate this for the case 
of the <j)-B cross-correlation; the il-E case will be similar. 
We define a new quantity C^ max , defined such that 

Ct B =AC* B nax . 
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If we assume that the instrumental noise on Cf 3 is Gaussian, so that is a random variable drawn 

from a Gaussian probability distribution with variance ^AC^ S ^ and mean A C^ ax , we can find the maximum- 



likelihood estimator for A to be 



s-l<p±5 

L ^ L L max 



(Aery 2 

El (c'lLx) 2 ( ac l B ) 2 ' 

Then, assuming that the instrumental noise is uncorrelated between different multipoles, the variance of this estimator 
is given by 



A' 



E (« ax ) 2 (ACf y 

. L 



The maximum signal-to-noise with which we can measure this amplitude is given by 



A v 



A 2 



E (« ax ) 2 (ACf )" 



1/2 



(52) 



The same method can be used to calculate the obtainable signal-to-noise from the Q-E cross correlation, giving 



Nj 



QE 



E 



jriQ,E \ 



1/2 



(53) 



Using the values of the cross-correlations and their errors calculated above, we find that the obtainable signal-to- 
noise from measurement of these cross-correlations is 0.002 for C^ E and 0.01 for Cf 3 . These numbers are too small for 
us to have any reasonable expectation of detection using the Planck satellite. Recalculating the above errors assuming 
an ideal CMB experiment, with no instrumental noise and infinite resolution, the values of the signal-to- noise only 
change by a factor of two, indicating that this method is not likely to be a promising way to detect a chiral GW 
background. 



CONCLUSIONS 



BiPoSHs are a formalism to describe correlations between two different spherical-harmonic coefficients of the CMB 
temperature field, which can occur if the CMB temperature field is not exactly Gaussian or statistically isotropic. 
This paper introduces odd-parity BiPoSHs, a set of BiPoSHs that has not yet been studied, and details how they can 
be estimated from knowledge of the CMB temperature fluctuations. 

We calculate the even- and odd-parity BiPoSHs that are sourced by gradient- and curl- type deflections of the CMB, 
respectively, and from this we obtain estimators for these deflections in terms of the BiPoSH coefficients. We show 
that lensing by scalar metric perturbations causes only gradient-type deflections, and thus only sources even-parity 
BiPoSHs. However, lensing by GWs produces both gradient- and curl-type deflections and thus sources both even- and 
odd-parity BiPoSHs. We calculate the expected power spectra of deflections due to scalar and tensor perturbations 
and their errors, and conclude that a reasonable signal-to-noise measurement of the amplitude of the GW background 
cannot be obtained from these autocorrelations even with the ideal CMB experiment, and thus from autocorrelations 
of the BiPoSH coefficients. 

Although lensing by GWs produces both even- and odd-parity BiPoSHs, their opposite parity implies that they 
could not be correlated. However, in the presence of parity-violating physics, such as a chiral GW background, this 
parity argument breaks down and we might expect a correlation. We consider such a cross-correlation, and encourage 
its measurement even though the likelihood of observing a cosmological signal is low. 

A GW background also produces signals in the E- and 5-type CMB polarization spectra, which are of even and 
odd parity, respectively. We consider the possibility that a chiral GW background would produce cross-correlations 
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between opposite-parity components of lensing and polarization, and calculate the expected magnitude and errors of 
such cross-correlations. Although we find that the likelihood of observing a cosmological signal is low, we encourage 
the measurement of these cross-correlations since such a detection would provide evidence of important systematic 
errors or even new parity-breaking physics. 

Here we have discussed BiPoSHs constructed from temperature multipole moments only, but the formalism can 
be generalized to include the polarization as well. It may also be that inclusion of the polarization improves the 
sensitivity to these parity-breaking, and other, signals. We plan to pursue this analysis in future work. 

Finally, we note that weak-lensing distortions of distant galaxies can also be decomposed into curl and gradient 
components j32l 153] . Similar tests for parity violation can thus also be carried out with weak lensing of galaxies. 
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